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ABSTRACT 

The formation of galaxies is regulated by a balance between the supply of gas and the 
rate at which it is ejected. Traditional explanations of gas ejection equate the energy 
required to escape the galaxy or host halo to an estimate for the energy yield from 
supernovae. This yield is usually assumed to be a constant fraction of the total available 
from the supernova, or is derived from the assumption of a consistent momentum 
yield. By applying these ideas in the context of a cold dark matter cosmogony, we 
derive a lst-order analytic connection between these working assumptions and the 
expected relationship between baryon content and galaxy circular velocity, and find 
that these quick predictions straddle recent observational estimates. To examine the 
premises behind these theories in more detail, we then explore their applicability to a 
set of gasdynamical simulations of idealised galaxies. We show that different premises 
dominate to differing degrees in the simulated outflow, depending on the mass of the 
system and the resolution with which it is simulated. Using this study to anticipate the 
emergent behaviour at arbitrarily high resolution, we motivate more comprehensive 
analytic model which allows for the range of velocities with which the gas may exit 
the system, and incorporates both momentum and energy-based constraints on the 
outflow. Using a trial exit velocity distribution, this is shown to be compatible with 
the observed baryon fractions in intermediate-mass systems, but implies that current 
estimates for low-mass systems can not be solely accounted for by supernova winds 
under commonly-held assumptions. 



1 INTRODUCTION 



Any viable theory of the formation and evolution of galaxies 
should be able to account for the mass of baryons contained, 
or rather not contained, in the massive collapsed regions 
that host galaxies. Observational constraints on the location 
of baryons in the Universe imply that the fraction within 
these 'halos' can be many times less than the cosmic baryon 
fraction, /j, « 0.17 (e.g. Komatsu 2011), and that the extent 
of the deficit is clearly dependent on the host's mass. This 
can be seen from the estimated baryonic and total masses 
from seven separate surveys which were collected together 
in one figure in the review by McGaugh et al. (2010); data 
which are reproduced here in our Fig. [I] 

The established explanation for this defecit, dating from 
long before such observational data were available, is that 
baryons can be driven from the galaxies - and their host 
halos - by supernovae explosions (Matthews & Baker 1971). 
This account is based on the premise that the energy re- 
quired to escape the galaxies' gravity is readily available 
from the supernovae. Because the gravitational potential 
barrier will increase with host halo mass, the fraction of 
the supernova-driven wind which escapes might intuitively 
be expected to be greater for lower mass systems, and this 



does indeed seem to be qualitatively upheld by the mass- 
dependance seen in the modern data. 

A more quantitative version of this theory was then de- 
veloped by Larson (1974), who equated this potential bar- 
rier with an estimate of the energy yield per supernova (and 
hence per mass of stars formed). In ij2] we review the ar- 
guments in this classic theory and, by updating the basic 
premises to include a cold dark matter component in the 
halos, show how it leads to the 1st order theoretical pre- 
dictions for baryon fractions which are overlaid with the 
observational estimates in Fig. [I] We also take the oppor- 
tunity to contrast the scaling expected from the traditional 
assumption, of consistent energy conversion to the ejected 



material (S 2.1 1, with the alternative working assumption of 



a consistent momentum yield (S2.2I 



We then go on, in [J3] to investigate how modern simu- 
lations of disk galaxies relate to these analytic theories, us- 
ing aspects of the theory to understand the behavior which 
emerges from these numerical experiments, depending on 
the parameters of the simulated supernovae and the resolu- 
tion with which the system is simulated. This investigation 
is then used to motivate revisions to the traditional premises 
and in fj4] we present a alternative derivation of the outflow 
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which is built on more realistic assumptions, whilst keeping 
to the spirit of the original elegant theories. 



2 THEORETICAL BACKGROUND 

2.1 Consistent energy conversion 

The earliest attempt to quantify the mass outflow due to 
supernova was by Larson (1974). This brought the idea of 
Matthews & Baker (1971) together with work on the evo- 
lution of supernova remnants (Cox 1972) to estimate that 
10% of the original supernova energy will be retained as 
thermal and kinetic energy in the g ELS, Et figure calculated by 
considering the energy that would still be contained within 
supernova remnants at the point at which they begin to 
overlap. 

This point naturally depends on the conditions in the 
interstellar medium, but by including two key properties 
(star formation rate and gas density) as variables in the 
derivation, the final estimate was shown to be sufficiently 
insensitive to variations in these conditions that it was con- 
sidered, "under most circumstances likely to be encountered 
in practice" , accurate "to within about a factor of 2" . 

This figure could then be incorporated into their model 
of gas ejection under the additional premise: 

"...that all of the available energy of 0.1-EgjM is expended in 
removing gas from the galaxy." (Larson 1974) 

This premise can be expressed as the assumption of con- 
sistent energy conversion, in that the energy yield is 
consistent from one system to another and is all converted 
to the gravitational potential energy gained by the gas in 
escaping the halo. This amounts to the relation: 



Host Halo Mass, M„ 
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(1) 



where M esc is the mass of ejected gas, v c is the characteris- 
tic velocity of the system (defined, for example, as the maxi- 
mum circular velocity) and n is the fraction of raw supernova 
energy, Esn, that is converted into the energy gained by the 
gas. Assuming the supernova yield to be proportional to the 
mass of stars formed, M*, this can be written: 
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The original calculations by Larson (1974) proceed from an 
effective energy release per supernova of 10 44 J and a mass of 
stars formed per supernova of IOOM0 , so the assumed yield 
of 10% corresponds, in our notation, to w§ N « 1O 41 J/M0 « 
(220 kms -1 ) 2 . 

The adopted figure for energy release of 10 44 J by Lar- 
son (1974) was attributed to observational estimates by 
Minkowski (1967), but has consistently appeared in the lit- 
erature since then without any citation at all. More recent 
justification for this figure has been provided by fitting evo- 
lutionary models to the observed variation in luminosity of 
the supernovae 1987A (e.g. Shigeyama & Nomoto 1990), 
1993J (e.g. Nomoto et al. 1993), and 19941 (e.g. Iwamoto 
et al. 1994), these having originated from stars with main 
sequence masses in the range 10-20Mq, and thus being rep- 
resentative of the majority of the supernova population. 

The scaling of eqn. [I] was subsequently adopted by 
White & Rees (1978) as a component in their analytic model 
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Figure 1. Galactic baryonic mass estimates from several publi- 
cations, shown as a function of the virialised host halo mass (es- 
timated from observed circular velocity or velocity dispersion). 
The lines show the simple prediction of self-similarity (dashed 
line), the prediction based on consistent conversion of supernovae 
energy to gravitational potential energy gain by the escaped gas 
(solid line and a similar prediction based on consistent momentum 
yield (dotted line) using the values v w = dsn = 300kms _1 . This 
plot is deliberately similar in layout to fig. 1 of McGaugh et al. 
(2010) and uses the same conversion from characteristic speed to 
host halo mass: M V /1O 12 M « (u c /187kms- 1 ) 3 . Round symbols 
represent rotationally supported disks, where v c is the outermost 
measurement of circular velocity. Square symbols represent pres- 
sure supported systems for which v c = %/3a is assumed, where a 
is the observed line-of-sight velocity dispersion. 



of galaxy formation, and the arguments behind it were re- 
visited and extended by Dekel & Silk (1986), though with 
rather detailed assumptions. White & Frenk (1991) felt these 
were difficult to justify, but they did adopt the scaling (if 
not the parameter value) in their models. Indeed, though 
different models employ their own particular definitions and 
contexts for its application, the assumption remains in use 
to this day (e.g. De Lucia et al. 2011). 

The mapping from supernovae physics to this observed 
trend in galaxy properties can be appreciated more directly 
if more approximately, by simply combining eqn. [I] with the 
additional constraint that the retained and escaped masses 
must add up to give the cosmic fraction, /b, of their host 
halo, (Mb + M esc = /bM v ). This leads to a simple relation- 
ship between remaining baryonic mass, Mb, and total virial 
mass, M v (or characteristic speed v c ) that is plottecQin Fig. 



m 



1 The line in Fig. ^ shows the locus of /* RJ 1. Real systems, 
where /*,$!, would therefore be expected to lie above this line. 
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So, placed in the context of CDM cosmogony, Larson (1974) 
seems to have predicted, to an order of magnitude, the stellar 
masses for a wide range of cosmic structures and, crucially, 
this prediction stems from a simple physical premise. 

But can basic energetic considerations really be a satis- 
factory explanation for this observed relationship? Can we 
really suppose that the converted fraction r\ (or equivalently 
wsn) is the same in such a wide variety of environments? 
Even if so, why would we assume further that this remaining 
energy, carried by the out-flowing gas, is then almost exactly 
exhausted in overcoming the gravitational potential? 



2.2 Consistent momentum yield 

The discussions reviewed so far have been focussed on the 
energy provided by supernovae. Yet, over time, stellar fu- 
sion generates a comparable amount of energy to the to- 
tal produced by these more powerful, but rare events. The 
reason that supernovae are held responsible for the mass 
loss is because the energy released in this process couples 
effectively to baryon^] With some fraction of the energy 
carried by massive particles, moving at some non-relativistc 
mean speed, v w , the momentum conveyed to the surround- 
ing baryons per mass of stars formed (psn ~ Esn/v w ) can 
be much greater than (or at least comparable to, see Murray 
et al. 2005) the total conveyed by the stars over their entire 
lifetime (p* £ E+/c where c >> w w ). So it is the momentum 
yield that is crucial for driving a galactic wind. 

An estimate of outward momentum may therefore seem 
a more relevant starting point for the derivation of the out- 
flow mass. This is further motivated because this quantity 
is understood to be roughly conserved (McKee & Ostriker 
1977) after the initial adiabatic expansion phase. So the 
need to continue modeling the remnant through the cooling 
phase, to derive a final energy conversion, is eliminated from 
this alternative lst-order estimate of the resulting outflow. 

Such an approach has some common ground with the 
one reviewed in §2.1[ as both contain the same assumption 
of energy conservatiorj^] between some known phase in the 
galaxy, and some point far from the galaxy where all the 
energy is supposed to have been expended in overcoming 
gravity. In both cases the working assumption is that, in 
order to have escaped, the gas of mass, M out , must at some 
earlier point have been moving outward with velocity v ^ v c . 

The difference in this second case is to associate this 
with a momentum, of about M out i>c, and suppose that this 
must (neglecting gravitational influence at this stage in the 
remnant's evolution) be less than or similar to the momen- 
tum yield at the end of the adiabatic phase. If it is this 
yield, per mass of stars formed, that is reasonably consis- 
tent across a range of systems (psn ~ M*« w , where u w is a 
mean system-averaged value for the specific momentum of 
the gas), this would imply: 



2 The key comparison being between the opacity at the centre of 
a collapsing massive star, and that in a stellar atmosphere. 

3 Both make use of an "energy conservation" argument, hence 
the need to take care when choosing terminology. 



The inequality also reflects the fact that this estimate of the 
ejected mass corresponds to its outward velocity tending to 
zero as it leaves the system (an approximation which will 
be redressed in sQ. As can be seen in Fig. [I] this reveals a 
much tighter constraint on the low-mass outflow that was 
apparent from the energy budget. However, it is important 
to remember that in the extreme low mass limit there are no 
stars at all, and hence no supernovae. So perhaps we should 
not expect this asymptotic behavior to be fully accounted 
for by supernova winds, and should recognise the additional 
importance of natural cooling thresholds and gas stripping 
(e.g. Stringer, Cole & Frenk 2010). 

Though derivations differ considerably, this alternative 
scaling has become established as being associated with the 
momentum-conserving properties of winds (we emphasise 
again that this conservation is only one aspect of the argu- 
ment) and has been applied when implementing supernova- 
driven winds (and indeed AGN-driven wind) in many pub- 
lished galaxy evolution models (e.g. Oppenheimer & Dave 
2006). 



2.3 Application to phenomenological models 

So far we have briefly reviewed two approaches which link 
escaped baryonic mass to the characteristic speed of the host 
system: the assumption of a consistent energy conversion of 
supernova energy to the outflow (eqn.[T| and the assumption 
that the momentum yield is consistent (eqn. [4| . A third ap- 
proach is to take a more general formula for the relationship, 
and tune the parameters to match simulations of the outflow 
process, and/or to match observations in conjunction with 
a galaxy formation model. 

Such an approach was taken by Cole et al. (1994), who 
analysed the earliest simulations of galaxies by Navarro & 
White (1993), where the effect of supernovae was explored 
by assigning velocity perturbations to the gas around star 
formation sites in their models. This implementation did not 
lead to the relationship given by ([I]), but the results could be 
matched by the more general expression (in our notation): 



Af out 



(5) 



where M out is the mass of gas which has been ejected from 
the disk (as opposed to the gas which has escaped the host 
halo, M csc , used in eqn. [I]) This empirical approach is not 
generally physically motivated (though of course it reduces 
in the cases a = 2 and 1 to the two approximations described 
above) . 

The free parameter a was found to increase monoton- 
ically with the fraction of energy assigned to the velocity 
perturbations in the simulations (up to a — 5.5 in the case 
of a 20% allocation). This general relationship between mass 
outflow and circular (or escape) speed was then adopted in 
many models, often fitting a by empirical constraints, and 
it remains in use (for the latest discussion, see Bower et al. 
submitted) . 

So, in the language of eqn. [5] a = 1, a — 2, and a whole 
range of other empirically motivated values are all being si- 
multaneously called upon to account for various different 
observational features of the galaxy population. They can 
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not all be correct at once. But arguing for any one simple 
approximation over the other is surely less progressive than 
seeking a more inclusive model which may reduce to certain 
simple scalings in particular regimes. Certainly, some align- 
ment of these apparently conflicting working assumptions is 
essential if any kind of consensus on this aspect of galaxy 
formation is to emerge for the literature. 

The first exercise towards this goal, presented in Sj3] is 
to update the work of Navarro & White (1993) and Cole et 
al. (1994) by investigating to what extent momentum and 
energy conversion are upheld by the results of the latest SPH 
simulations of the process. This study is integrated with an 
investigation into the modelling technique itself, consider- 
ing variation of both wind parameters and the numbers of 
particles used to represent the systems. Guided by the impli- 
cations of these numerical experiments, we then present in 
Sj4]a revised analytic picture, which reflects the reality that 
gas will leave its host galaxy with a range of exit speeds, and 
includes both energy-conserving and momentum conserving 
constraints. 



3 FEEDBACK IN IDEALISED DISK 
GALAXIES 

3.1 Summary of the numerical experiments 

The simulations of idealised, isolated disk galaxies that will 
be the focus of this section are generated from the same 
initial conditions as used by Dalla Vecchia & Schaye (2008) 
(after Springel, Di Matteo & Hernquist 2005). We follow 
their evolutiorj^Jfor 0.5Gyr using the modified version of the 
GADGET-3 code (last described in Springel 2005) . 

This version has been used in the Overwhelmingly 
Large Simulations (OWLS) project (Schaye et al. 2010) and 
in the Galaxies-Intergalactic Medium Interaction Calcula- 
tion (GIMIC, Crain et al. 2009). In particular we use the 
'reference model' of OWLS. Briefly, sub-grid physics is im- 
plemented such that stars form at a rate which guarantees a 
Kennicut-Schmidt law for a disc in hydrostatic equilibrium 
(Schaye & Dalla Vecchia 2008) and gas cooling includes the 
contribution from hydrogen, helium, 9 'heavy elements' and 
bremsstrahlung, as described in (Wiersma et al. 2009). We 
assume the gas to have a solar abundance pattern with so- 
lar metallicity, and stars are assumed to form with the IMF 
proposed by Chabrier (2003). Cosmological simulations us- 
ing this code form plenty of galaxies that are similar to the 
Milky Way, in terms of stellar (Font et al. 2011) and X-ray 
properties (Crain et al. 2010). 

The idealised galaxies have dynamic baryonic compo- 
nents, with mass fraction 0.054, embedded in static dark 
matter halos containing the remainder of the total mass, 
Mhaio, with mass distribution M(r) = Mr /(r + raq) 2 
(Hernquist 1990). The lengthscale, ruq, for each galaxy is 
given in Table |3.1| The baryonic mass is divided between 
a stellar bulge with total mass fraction of 0.014 (also dis- 
tributed with a Hernquist profile) and a disk of fraction 0.04. 

4 This time interval was considered to be appropriate for the 
purposes of this study, being long enough to span many dynami- 
cal times for all the systems, whilst being short enough that the 
isolation from the cosmology was an acceptable approximation. 



Table 1. Parameters of the idealised galaxy simulations. 



Total Mass/M 


10 10 


10 10.5 
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10 11.5 
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Galaxy Mass/M Q 


5.4- 10 s 


1.7- 10 9 


5.4 ■ 10 9 


1.7- 10 10 


5.4- 10 10 


Disk ro/ * pC 
ZD/kpc 


0.53 
0.053 


0.78 
0.078 


1.1 
0.11 


1.7 

0.17 


2.5 
0.25 


Bulge r Hq /kpc 


0.053 


0.078 


0.11 


0.17 


0.25 


Halo rHq/kpc 


6.52 


9.56 


14.0 


20.6 


30.3 



Of the disk mass, 70% is contained in a pre-existing 
disk of old star^] and the remaining 30% is composed of 
cold disk gas. Both disk components have surface density 
distributions which decrease exponentially with radius with 
the same scalelength, set iteratively such that the angular 
momentum of the disk is a particular fraction of the total 
(with spin parameter, A = 0.033, see Springel 2005). The 
disk scaleheight and bulge scalelenghts are then both set to 
a tenth of this value, with the stellar disc following the profile 
of an isothermal sheet and the vertical gas distribution set 
up in hydrostatic equilibrium using an iterative procedure. 
The simulations were run with a default gravitational soft- 
ening length of lOpc, chosen so that this is always smaller 
than the Jeans length at the highest surface densities we 
expect]^] 

The effects of supernovae are simulated by assigning 
particles surrounding a star forming particle a velocity kick, 
of magnitude v v , in a random direction. Neighbouring par- 
ticles are chosen at random such that their total mass is 
proportional to the initial mass of stars formed. The default 
choice for this "mass loading" in these models is: 

AM W = 2AM* , (6) 

where A denotes a change over a given time interval. To 
yield an appropriate total energy, this mass loading cor- 
responds to an initial wind kick, « w = 600 km s -1 (giving 
AEsn « 7 x 10 41 (AM*/M Q )J). These fiducial values follow 
Dalla Vecchia & Schaye (2008), and are adopted here when- 
ever other parameters in the model are under scrutiny. It 
is important to note that all the available supernova energy 
is allocated to these velocity perturbations, the radiation 
losses being left to emerge from the SPH calculation itself 
(rather than assumed upfront). Analysis of this emergent 
energy conversion is the focus of the following section, |3.2[ 
which begins our analysis of these numerical experiments. 

We note that this is by no means the only possible 
numerical experiment that can guide our understanding of 
these processes. By simulating only a small section of the 
interstellar medium, it is possible to resolve length scales 
which, in the implementation here, have to be left to the 
sub-grid configuration of the model. Such investigations are 
in progress (?) but these smaller volumes must necessarily 
make other assumptions regarding their interaction with the 

5 In this model, the contribution to the total energy yield from 
ongoing winds or supernovae from this "old" population is as- 
sumed to be negligablc. 

6 Note that since we integrate over a short time in a galaxy with 
a significant dark halo, two-body relaxation effects are mostly 
negligible. 
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system in which they are supposed to be imbedded. For this 
study, we choose to focus on simulations of entire galactic 
systems. 

3.2 Initial comparison with consistent energy 
conversion argumnets 

To connect with the review of analytic treatment in EjlJ 
we begin by presenting a few pertinent results from five 
idealised disk galaxies covering a range of host halo (and 
galaxy) mass. The results are plotted in Fig. [2] as a function 
of characteristic velocity, v c , here defined as the root-mean 
absolute specific energy, v c = y/—eo, measured for all the 
gas in the galaxy (after several dynamical times, but be- 
fore perturbation by supernovae). This working definition 
is adopted as it is a well-defined quantity in the simulation 
and is consistent with its application in the analytic model 
(which is concerned with the energy deficit that has to be 
overcome). Also, this and the observable v c are mutually 
satisfactory proxies. 

The top panel shows the mass, M CBC , which has "es- 
caped" the galaxy (outside the disk and with positive total 
energy). The lower panel shows the totaQ energy of this 
body of gas, the solid points being the value at the final 
time-step. To provide a context for these "excess" ener- 
gies, a faint line is included showing the approximate scaling 
that would result from a common residual wind velocitjrl 
= 140 km s -1 , and the mass scaling of equation |TJ. 

For now, the most important thing to note is that these 
measured residuals of final energy are just as significant, if 
not more so, than the initial energy of the same body of 
gas, shown as open points in the lower panel. So in this par- 
ticular model, at least, it seems that there is much more 
energy expended than is required to just "remove gas from 
the galaxy" , as described in the quotation in §2.1| A signifi- 
cant fraction is expended in the form of kinetic and thermal 
energy of the gas as it escapes. In the case of less massive 
galaxies (where the wind velocity turns out to be several 
times greater than the characteristic velocity) it is actually 
the gravitational energy that could be neglected, if anything. 

This can be clarified by reference to an extended version 
of |TJ which identifies some other energy sinks: 

1 



-Esn = E\o 



+ ^M esc v; 



+ E 



trapped 



+ E E 



(7) 



Equation |TJ implicitly assumes that the final two terms here 
can be neglected. As just discussed, these particular results 
clearly lend no support for neglecting the excess energy, E ex . 
Indeed, in reality, it is difficult to see how the mass of gas 
affected by the supernovae could be naturally tuned such 
that it will always coast gently out to rest in the intergalactic 
medium. 

Conversely, it also seems probable that plenty of gas 
will acquire some energy from supernovae but not enough 
to escape (corresponding in eqn. rapped > 0). In order 

to satisfy energy balance of the system over many dynamical 
times, there is an argument that this must be eventually lost 

7 Kinetic, potential and thermal energy are all included. 

8 Direct measurements of this reveal true mean values in the 
range D ex ~ 100 to 150kms _1 , increasing with v c . Exit velocities 
are investigated in 93.41 
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Figure 2. A comparison between traditional feedback theory and 
results from state-of-the-art SPH simulations. The points in the 
upper panel show the escaped gas mass, M eBC , in five idealised 
disk galaxies after 0.5Gyr of evolution (in units of the the mass 
of new stars formed, AM*), as a function of their characteristic 
speed, v c . The galaxy masses range from 5 x 10 8 to 5 x 10 10 M Q . 
The lower panel reveals the energetics. Solid points show the to- 
tal energy that is escaping with the gas and open points show the 
initial energy of the escaped gas (after a few dynamical times, but 
before perturbation by supernovae). Both are plotted as a frac- 
tion of the total energy input. The dotted line in both panels 
attempts to apply equation to these results. In the bottom 
panel shows the mean of the measured efficiencies, n. The top 
panel puts this value back into |T]| to see if it agrees with the 
results. The faint solid line in the lower panel shows the approx- 
imate outgoing energy which would result from the mass scaling 
in this equation, together with a common mean exit speed 4 , v ex - 



(and hence covered by E\ ost ). But, even if so, this "fountain" 
process must still be considered in any reasonable derivation 
of the total energy loss. 

Yet, despite all these caveats, the measured mass out- 
flow in the model, shown in the upper panel of Fig. [2j does 
appear to agree with |T]). The reason for this can be quickly 
appreciated: Though the total energy fraction converted is 
variable (as discussed above) , the potential energy which has 
been overcome to escape is similar in all systems. So the 
assumption that it is the only energy sink, though in fact 
incorrect, appears to work. 

This is shown by the open points in Fig. [2] together 
with a dotted line indicating the mean value, n. When this 
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valu^jis used in Jy, the small scatter around it means that 
there is reasonable agreement with the measured mass in 
the upper panel. 

So here we have a plausible scenario in which a fre- 
quently adopted scaling does, broadly, apply, but in which 
the usual assumptions are found to be only partially true. 
This is an illustration of the dangers of confirmation by af- 
firmation. However, an illustration is all that it can be at 
this level of analysis. The degree to which any particular 
implementation of winds, or set of numerical methods may 
represent reality is far from established. What we can try to 
establish is this: 

• The physical cause for certain aspects of the model's 
behaviour. 

• Which aspects therefore seem unreliable? 

• Which aspects are indeed telling us about real galaxies? 

• How can our physical understanding of the process be 
modified in the light of this investigation? 



3.3 Variation of injected wind momentum 

To begin analysing the modelling method, we briefly present 
the effect on the outflow of varying the wind speed and mass 
loading. To help isolate the effect, the total wind energy in- 
put, -Esn = §Afwfw! is kept constant for all the realisations 
in this test. 

Fig. [3] shows the mass, energy and momentum out of 
the plane that escape the same galaxy under five different 
initial wind speeds, but with the same total energy budget. 
With high wind speed, to the left of the plot, a small number 
of gas particles (about ten thousand per timestep) receive a 
large velocity kick. At these high speeds, the change in speed 
due to gravity is small and almost all of the momentum out 
of the plane is conserved in the outflow. 

With low wind speed, towards the right, a larger number 
of particles (a few hundred per timestep) receive a small 
kick, and almost all the energy and momentum is absorbed. 
However, both the energy and, crucially, the mass which 
escapes are quite robust to changes in mass loading at the 
highest values of this parameter. It should be noted that the 
default choice (« w = 600 km s" 1 ) is in the range where the 
mass outflow is sensitive to this value, for this halo mass. 



3.4 Convergence and limiting behaviour 

In addition to the importance of parameter choice, galaxy 
simulations such as these can also be very sensitive to res- 
olution. In this section, we illustrate the differing effects of 
the same physical assumptions applied to initially identical 
model galaxies which are resolved at a range of different 
values for gas particle mass, m p . 

For a clear study of the effects of varying particle mass, 
the same gravitational softening length was used for all the 
runs. Different values would have been adopted for different 
resolutions if the simulations were being run in their own 
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Figure 3. The effect of assuming different intrinsic mass load- 
ing in the modelling of supernova-generated wind, given a fixed 
energy yield for a halo of Mh a i = 10 12 Mq. As the increased 
mass loading leads to a higher initial momentum, this quantity is 
plotted as the main x-axis (the wind speed decreases from left to 
right, as shown at the top of the figure). The top two panels show 
the same quantity as those of Fig[2] The bottom panel shows the 
momentum of the gas out of the plane of the galaxy, as a fraction 
of the initial momentum in these directions. 



right, rather than as part of a convergence test. A differ- 
ent choice would have had a lesser effect than varying the 
particle mass (see Appendix |A| . 

For this more detailed investigation, we adopt a slightly 
more precise definition of ejected gas than was used in [JT] We 
define two planes which initially enclose all the particles at 
the highest resolutiorp°]by a margin of 10%. (So the outflow 
boundary exceeds the unperturbed height of any particle by 
a factor of at least 1.1). It turns out that this is almost the 
same for all systems, between 2 and 2.2kpc above and below 
the disk plane. A number of useful quantities can then be 
measured from each simulation, for any chosen time period 
up to the total duration of the calculation: 



10 The highest resolution was limited by the number of particles 
9 In this case n = 2.4% and -Esn = (600 km s^ 1 ) 2 M t , so = (rather than particle mass), and a practical maximum number 

93kms — 1 . turned out to be, for gas particles, N p = 2.3 X 10 5 . 
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• The gas mass which flows out through these two bound- 
aries, Mout- 

• The distribution of this mass as a function of velocity 
component perpendicular to these planes, v z . 

• The mass of gas which flows back through the boundary, 
having been outside it ("re-cooled") 

For all of these components, it is instructive to distin- 
guish between "wind particles", that were directly assigned 
kinetic energy as a result of being near a simulated su- 
pernova, and particles that have been indirectly affected. 
The total mass of directly effected wind particles is denoted 
A^fwind- From the geometry of the randomly orientated kicks 
that they receive in this implimentation, the initial distri- 
bution of velocities out of the plane is dA//dju z j = M w /ij w 
and the initial momentum away from and perpendicular to 
the plane is |M w v w . 



3.4-1 Dwarf galaxy 

Firstly we focus on a dwarf galaxy with Afhaio = 1O 1O M0, 
where the effects of using different resolutions are particu- 
larly palpable. Images of the gas in the simulation at two 
different resolutions are shown in Fig. [4] and some quantita- 
tive results from simulations of this system, across a range of 
resolutions, are shown in Fig. [6] The top panel of this figure 
shows the total gas mass which has been directly accelerated 
by simulated supernovae (the "initial wind"), the total mass 
which has been ejected beyond the nominal galaxy bound- 
ary (defined above) and also the re-cooled mass which has 
returned to the inner region after being ejected. 

These cumulative values, over 0.5Gyr of evolution, pro- 
vide a simple illustration of the behaviour of the simulation, 
which changes markedly as more and more particles are used 
to represent the galaxy. At the very lowest resolution, all the 
supernovae in the galaxy during a timestep are represented 
by a single gas particle which is allocated all the energy and 
duly exits the disk. So in this limit, M out — > M v - ln d. 

This can be appreciated by reference to the middle row 
in Fig. [6j which shows the time evolution of three of the 
cumulative quantities from the top panel. In the leftmost 
panel, the ejected mass fraction can be seen to be increasing 
in steps, corresponding to single wind particles exiting the 
disk regiorp*"! Gas which is not directly kicked by supernovae 
acquires very little of the energy and does not escape the 
disk. 

With more particles, and more interactions, such sec- 
ondary outflow does become possible. Because of this, the 
evolution of the system looks very different to the right of 
Fig. [6] compared to the left. As the number of particles in- 
creases (left to right), the initial supernova momentum is 
shared by a greater gas mass. So the total energy which the 
outflow carries out is less, but the mass loss is greater. To 
accompany this explanation, a faint line is included in the 



11 The difference between ejected and wind mass in this case 
is due to a single particle which acquired an energy kick at t ~ 
O.lGyr but remained contained in the disk. All the others escaped, 
so the two lines differ by the mass of one particle. The one particle 
happened to be contained because of its initial trajectory towards 
the disk centre. 



upper panel of Fig.[6]which shows the momentum conserving 
case: 



(8) 



where is the mass-weighted mean of the exit velocities. 

At low and intermediate resolution, interactions with 
the remaining disk gas are not generally sufficient to bring 
gravitational deceleration into play (i.e. v\ » v\ for much 
of the ejected gas) so equation Q, which neglects gravita- 
tional influence, is a good fit to the results in this regime. 
As the resolution is increased and exit speeds are reduced, 
more and more mass is ejected from the system. Eventually, 
at the highest resolutions, the speeds are reduced to v z ^ v c , 
gravity has a significant influence, and the results begin to 
diverge from the simple prediction of |8J. 

So, the broad behaviour of this numerical experiment 
can be understood in terms of physical arguments, and the 
principal results - star formation and outflow - do appear 
to be reasonably robust to changes in resolution above some 
reasonable particle mass (~ 10 4 Mq). This is supported by 
the behaviour shown in the top panel of Fig. [6] However, 
though the results may be consistent over this high-end 
range of resolution, the changing velocity distribution in the 
bottom row indicates that the calculation has not fully con- 
verged on all scales. 



3.4-2 Milky Way-like galaxy 

We now turn to the case of a larger galaxy, with Mi la i = 
1O 12 M0. Images of two realisations of the system can be 
seen in Fig. [5] and the key results are shown in Fig. [7] (in 
the same format as Fig. [6j| . 

To understand the behaviour of this larger system, we 
begin with the low resolution limit where interactions be- 
tween wind particles and the rest of the disk gas are negli- 
gible. The outcome is nearly as straightforward as the low 
resolution limit in the dwarf galaxy, where almost all wind 
particles simply escape. But because this larger system has 
at least some gravitational influence on the wind particles 
(v c ~ v w ) a significant fraction are retained even in the ab- 
sence of significant additional interaction. This statement 
implies something along the lines of: 



Mou 



AM 



dv z 



(9) 



For the simple low-resolution case of a flat initial mass distri- 
bution in v z , which results from randomly directed velocity 
kicks, equation (J9j) becomes: 



M 



M^ 



Vent 



for 



AM 

dv z 



M a 



(10) 



A dot-dashed line is included in Fig. [7] to illustrate this re- 
lation, using Vcut ~ Vc as a nominal value for the velocity 
required, on average, for gas to escape. 

Of course, the distribution of wind velocities is only flat 
at the lowest resolutions, so this does not match the results 
of the high resolution runs. As the system is modelled with 
more and more particles, the velocity distribution tends to 
the smooth form that can be seen in the bottom right panel. 
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Figure 4. Two realisations of the 1O 1O M0 halo, showing the projected gas density (brighter regions show higher densities). The left 
panel shows a low resolution version (which starts with 560 gas particles), corresponding to the second from left column in Fig . |6| The 
right panel is the highest resolution that was achievable (750,000 gas particles), corresponding to the far right column in Fig. |6] Both 
show the system at the final timestep, after 0.5Gyr of evolution. 




mp = 1.3 x 10 4 M o 



8kpc 



Figure 5. Two realisations of the 1O 12 M0 halo, in the same format as Fig. [4] but at a larger scale (as shown). The selected resolutions 
correspond to the second column (93,000 gas particles) and fourth column (3 million particles) in Fig. [7] In this figure, and in Fig|4] each 
particle is represented by a smooth radial density profile which reduces to zero at a distance of 2s, where s is the lesser of a stipulated 



maximum, s r . 



and the distance to the 16th nearest neighbour (s n 



[1.5,6.5]kpc for M halo = [10 10 , 10 12 ]M Q ). 
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Figure 6. The effect of resolution on the mass outflow and evolution of an idealised, isolated dwarf galaxy (halo mass= 1O 1O M0, baryonic 
fraction=0.054). The same initial conditions have been evolved using the GADGET code but with 15 different choices of particle mass. 
The gravitational softening length in all cases is lOpc. The middle panels show the time evolution of four of these, showing the gas 
mass fraction which has escaped (solid line), and the fraction which has escaped from the disk but returned (dashed line). Also shown is 
the fraction which has been directly allocated wind energy (dotted line), which in this set of simulations is just M w = 2M*. The upper 
panel shows the final value of each of these quantities as a function of the particle mass used in the simulation run. Added to this panel 
is a faint dot-dashed line showing the result of straightforward conservation of initial momentum out of the plane |8]) . The lower panel 
shows the distribution of escaped gas as a function of the outward velocity component, v z , at the point when it crosses the nominal 
galaxy "boundary" (see main text). The vertical dot-dashed line shows the mass-weighted mean, (in kms -1 ) which is used in cqn.|s] 
Note that the integral under the total ejected mass is equal to M ou t/M w ; n( j, but the shaded area (which includes only wind particles) is 
not equal to unity because some wind particles do not escape. 



10 M. J. Stringer et al. 



Increasing resolution > 



0.8 
o 0.6 

m 

do 

^ 0.4 
0.2 



1 1 1 1 1 1 
halo © 


1 1 1 


| 1 1 1 1 | 1 1 1 1 | 1 1 1 










— , 


a ..-a Inital wind 

"" E1 --.. ..a-"' ~ 
■ a- " 












Ejected 


i i i i 1 i 


i i i 


1 i i i i 1 i i i i 1 i i i 



o ■— ^ 

8 i 



7 



^ l 
o 

^ 0.8 

EG 
CO 

0.6 
^ 0.4 
^ 0.2 

3 
O 

S 



V 



log( Particle Mass /'M ) 

V V 



I 4 
l 

V 



I M p =5360.4 xlO 4 M Q : 


M p =42.9 xlO 4 M Q I 


M p =10.7 xlO 4 M Q : 


M p =1.3 xlO 4 M : 


.j-^i i i 1 — >■ —r T i 


.x^T, , _L -r T r- 


.!>< i j _4 -T "T i i 


- 1 — \ ~t i "|" i""""* 



0.25 



0.25 0.25 

Time / Gyr 



0.25 




Figure 7. The effect of resolution on the mass outflow and evolution of an idealised, isolated Milky Way-like galaxy (halo mass= 10 12 Mq, 
baryonic fraction=0.054). The key is the same as Fig.|6] except that the faint dot-dashed line for the limiting behaviour is different, as 
explained in j]3.4.2| Note, again, that the mass of stars formed is just half the mass of the initial wind (dotted line). 
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Notably, at the very highest resolution, the distributions of 
directly and indirectly affected particles are beginning to 
look quite similar. This is an encouraging sign that the cal- 
culation is beginning to converge, but the total mass outflow 
is still lower with each successive increase in particle num- 
ber. 

Concern over the lack of convergence in Figs. [6] and 
[7] could be compounded by concerns for the particular im- 
plementation of winds that has been used in the models. 
Idealistically, one could argue that the precise choice of im- 
plementation would not matter with a sufficiently reliable 
SPH calculation, ejecting energy as random velocity per- 
turbations to the particles being in principle equivalent to 
updating the "thermal" energy of each particle by the ap- 
propriate value. But, in practice, such concordance is not 
the case without a more efficient time-integration scheme, 
as shown by Durier & Dalla Vecchia (20f2). 

Some simple tests, such as shown in Appendix [B] sug- 
gest that the consequences for these particular experiments 
(even in the case of the low resolution runs) is fortunately 
not critical. But ultimately, this is yet another reason to 
remember that the results presented in this section are an 
investigation rather than a test. That said, we return to our 
premise at the end of §3.2| If we can at least understand the 
behavior of our models (incomplete though they may be) 
in a physical way, we can use this to help motivate a more 
comprehensive understanding of real systems (even though 
the model itself may not perfectly represent them). 

Thus, if |9| works well for the simple velocity distri- 
butions seen in Fig. [7] at low resolution, it might well be 
expected still to apply at the highest available resolution 
and indeed in reality, if only we knew the appropriate veloc- 
ity distribution to use. So, if we could anticipate the velocity 
distribution towards which the simulated winds are begin- 
ning to converge, we would equally well be able to predict 
- and understand - the mass outflow that would emerge at 
arbitrarily high resolution. This prospect will be examined 
infl] 

3.5 Mass dependence 

The goal set up by iJTJwas to acquire a better understand- 
ing of the correlation between outflow efficiency and host 
galaxy properties, in particular the mass (or circular veloc- 
ity). To reinitiate this discussion, Fig. [8] presents some key 
results from the idealised simulations as a function of host 
halo mass, but now showing not only the highest available 
resolution that was attained for all systems (using the same 
number of particles) , but also the equivalent results at very 
low resolution. Reassuringly, the star formation appears to 
be robust to variation in resolution, so this discussion can 
focus on the behaviour of the outflow 

Turning to the outflow mass, the lower resolution case is 
easily understood by noting that, at all halo masses, the out- 
flow mass is approximately double the mass of stars formed 
(clear from the upper panel) , which is simply the "mass load- 
ing" employed in this feedback prescription (JsJ) . Broadly, all 
kicked particles escape, others are unaffected. 

At the highest possible resolution, the mass outflow is 
highly dependent on the local environment. As discussed in 
§3.4| this variation is the result of two principal competing 
effects: 
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Figure 8. The outflow of gas from idealised galaxies of different 
total mass, Afjialo' The top panel plots the ratio between the 
ejected mass and the mass of stars formed. The middle panel 
shows the fractions which have been ejected from the system and 
re-cooled and the lower panel shows the fraction of the initial 
gas mass that has formed stars after being evolved for 0.5Gyr 
(circles). In all panels, the highest resolution is shown as solid 
lines/points and the lowest choice as dashed lines and open points. 

• the distribution of energy across a greater gas mass 
(tending to increase the ejected mass due to conservation of 
momentum) and 

• the radiation of energy that such distribution will pro- 
duce 

(tending to reduce the ejected mass; lower velocities are in- 
sufficient to escape gravity). 

A quantitative account of the relative importance of these 
two influences will be the focus in Q 



4 TOWARDS AN IMPROVED PICTURE 

Armed with a better understanding of the simulations, we 
return to the task motivated in §3.2| to seek a link between 
the external effect of supernova feedback (e.g. M ou t), the 
internal conditions in the galaxy (e.g. vt) and the governing 
parameters (e.g. u w ). Some steps towards this have already 
been taken in limiting cases. 

The case of a neglible potential barrier, but significant 
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energy losses, was satisfactorily accounted for by a momen- 
tum conserving argument (JsJ) . This was compatible with the 
outflow in the limiting case (y%, » v%) a good approxima- 
tion for the lowest mass galaxy up to quite high resolution 
(Fig. |6| and also for the most massive galaxy when simu- 
lated with a very high choice of intrinsic wind speed (Fig. 

The case of a significant potential barrier, but negli- 
gible energy losses, was successfully accounted for by Q. 
When applied to the limiting case of ballistic wind parti- 



FT 



cles, perturbed only by gravity, (10 1, this accounted for the 
low- resolution limit in the high mass galaxy (Fig. [7]). 

The two ideas can be combined. To do this, we con- 
sider a momentum-conserving phase, leading to some dis- 
tribution of outward velocities AM/Av z , of which only the 
high-velocity tail, v > v cut is able to escape. This implies: 

nj P SN f°° AM a fill 

M out = / dv z (11) 

v z / dv z 

•> "cut 

where psn is the initial momentum provided by the super- 
novae in question. A more careful derivation considering all 
directions, rather than our simplistic consideration of out- 
flow directly out of the plane, would be of similar form but 
contain structural parameters. 

Further illustration of this approach can be provided 
by adopting a nominal form for the velocity distribution, 
dM/dv z , which appears in (111. We consider the simple 
form: 



dM 

dv z 



(12) 



where \M m v w has been substituted for the initial momen- 
tum provided by the supernovae. 

This ansatz for the converged velocity distribution can 
be compared with measurements from the highest resolution 
simulations that are available from our range. This compar- 
ison is shown in Fig. [9] which plots the velocity distribution 
of the ejected gas from three systems of different mass, sim- 
ulated with the highest number of particles with which all 
three could be completed. 

The v z that are measured from the simulations are 
the eventual exit velocities. Because ( 12 1 is an attempt to 



characterise the distribution of velocities at an intermedi- 
ate stage, these values are adjusted by the potential gain, 
A<J>j, at the point of measurement from the pre-kicked po- 
sition: ii at jj = Vv'i + 2A$i. This is somewhat crude but is 
more faithful to the concept behind eqn.[l2]and is adequate 
for the main purpose of this figure, which is to illustrate its 
application. 

The distribution of outflow g function of these 

adjusted velocities shows that only the high velocity tail 
makes it out of the disk, and the length of that tail will 
depend on the potential barrier in question. We model this 
as a smooth distribution cut-off at at some value, v cu t , which 
is determined by the depth of potential well from which the 
gas has escaped (and thus expected to be close or equal to 
the characteristic velocity, v c ). The cut-off is not abrupt in 
reality, but the approximation of an abrupt cut off may be 
an effective modelling strategy. 

As shown in Fig. |9j v cut does indeed scale convincingly 
with halo mass and v c . So the gravitational potential barrier 
presented to the supernova wind, known to be important 
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Figure 9. A comparison between the approach to modelling out- 
flow described by equation and the results of the controlled 
simulations of disk galaxies, run with the same number of particles 
(fb 2 X 10 5 gas particles). The upper panels show distribution of 
ejected gas as a function of outward velocity, from the simulation 
results (solid line) and according to eqn. |12| (dashed line). The 
mean, v, of the velocity distribution is fitted to match the 
simulation results and the cut-off velocity, v cu t , chosen to enclose 
the appropriate gas mass. These cut-off velocities are then plotted 
in the bottom panel as a function of host halo mass and com- 
pared with v c . The dotted line shows the theoretically expected 
scaling with M^i a . 
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Figure 10. A comparison between the simulation results, at the 
highest achievable resolution, and eqn. |13| which attempts to ex- 
trapolate beyond the resolution limits. This choice of v is based 
on the velocity distribution which best fits the outflow from the 
1O 1O M halo, so this can be interpreted as a prediction of the re- 
sults that would be found if all the systems were simulated with 
the particle mass that is possible for this lowest mass galaxy. 



since the work of Matthews & Baker (1971), can still be 
included in a simple estimate of the ejected mass, but in a 
context which is more faithful to the physics of the process. 



4.1 Reconciliation with S[T] 

As illustrated in Fig. [9] eqn. [Tl] can yield an expression for 
the outflow in terms of two parameters, both of which have 
direct physical significance: 

• The distribution of speeds for supernova shock fronts 
at the end of their momentum-conserving phase. 

• A cut-off velocity, « cu t, which approximately selects 
from this distribution the gas mass which will escape the 
high-density region of a galaxy (expected to be proportional 
to v c , or correlate with total mass as 



Using the example velocity distribution in (121, this gives 
an alternative equation for the mass outflow: 



M out 


2 v w 


M w 


V 2 


M out 


v w 




a — e 

V 



-V z / V 



dv. 



t/v 



(13) 



This is compared in Fig. [To] to the results of the highest 
available resolution runs of the five galaxies (the same points 
as Fig. [sj using values of v = 40 and D cu t = v c . 

The chosen v is the value that produced good agreement 
with the 1O 1O M0 halo, which could be run with the highest 



resolution (lowest particle mass). Thus, the line in Fig. 10 
can be viewed as an attempt to extrapolate from this to 
predict the outflow that would be found if we could resolve 
all the systems to this degree. 

For this reason, it is not agreement with the simulation 
results that is sought from Fig. [To] The goal is to supersede 
them, and to promote a more digestible explanation for ob- 
servations such as the dependence of feedback efficiency on 
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Figure 11. Application of the revised theory to relevant 

data from Fig. ^ The solid line shows the equation given, using a 
mean wind velocity of v = 70kms — 1 and a mean specific initial 
wind momentum per mass of stars formed of -u w = 1400kms _1 . 
The dotted line shows the limiting maximum outflow for this 
particular wind velocity distribution and value of ti w . 



gravitational potential in real galaxies. With this in mind, 
we return to the issue presented at the start of the paper to 
see how our analysis of the simulated galaxies might apply 
in the context of real ones. To begin with, we identify again 
the two natural limits in the theory. 

The simplest of these is at high mass, where the escape 
speed of the system is simply too high and supernovae do not 
generate the required wind speeds. This is the same physical 
limit in both theories (eqn. [I] and eqn. 111. 

At the very lowest mass, the two models differ. In eqn. 
[l] the outflow velocity is implicitly set to be zero at large 
radii and the only limit on the outflow is from energy con- 
servation. But in order to escape at all, the wind must have 
had at least some outward momentum, of order M ou tV c , and 
this momentum must have come from somewhere, implying 
M ou t ~ Psn/vc - The model encapsulated by eqn.lllldoes not 
ignore this, and there is indeed an equivalent maximum out- 
flow for a given momentum yield, psn, from the supernovae. 
For the velocity distribution of eqn. (12} this occurs at v = v c . 

At intermediate masses, there will be a gradual tran- 
sition between these two limiting physical cases. This is il- 
lustrated in Fig. |11| which completes the paper by applying 
the alternative theoretical approach alongside the relevant 
observational data from Fig. [T] which originally motivated 
the investigation. As indicated in the labelling, the baryon 
content of the lowest mass structures will be severely limited 
by additional factors, principally natural cooling thresholds 
and the stripping of hot gas by larger neighbours. A detailed 
explanation of the way these effects function in the context 
of hierarchical formation can be found in Stringer, Cole & 
Frenk (2010). For this reason, the low mass observational 
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estimates are not included. The high-mass structures, their 
baryonic contents also understood to be limited by processes 
other than supernovae, are similarly omitted. 

The values of « w and v in Fig. [IT] are chosen to illus- 
trate that the theory is compatible with the observations 
using physically appropriate parameter values. The value 
for the momentum yield^Jof psn = 1400M* kms -1 is some- 
what higher than the nominal value plotted in Fig. [l] but 
the locus is similar. This is because the models reviewed in 
[Jl] consider only one, characteristic value of momentum or 
energy for the wind, and assume that all the wind escapes. 
The revised theory allows for a distribution of velocities in 
the wind; the slow-moving fraction being retained by grav- 
ity, and the high-velocity tail escaping with, in general, a 
significant asymptotic kinetic energy. We consider this to be 
a more realistic picture. 



5 SUMMARY 

In this paper, we began by reviewing the 40 year old theory 
of supernova-driven outflow from galaxies, distinguishing be- 
tween the premise of a consistent energy yield and that of a 
consistent momentum yield. If adapted to CDM cosmology, 
we demonstrated that such theories can be used to derive 
lst-order predictions of the expected baryon content of cos- 
mic structures which are in agreement with observational 
estimates spanning many orders of magnitude in host halo 
mass (Fig. [TJ. 

With the aid of modern computational resources, we 
then scrutinised the arguments in these theories, following 
the extent to which certain premises apply as parameter 
values (Fig. [3f and particle number (Figs. [6] and [7| were 
varied in the simulation, identifying regimes dominated by 
momentum conservation and those where energetic limits 
take over. One traditional premise which was not upheld in 
these models was the supposition that all the energy con- 
veyed to the gas would be exhausted in overcoming gravity. 
The simulated outflows left the system with a range of veloc- 
ities, comprising a fraction of the total energy comparable 
to that expended against gravity. 

The tendency of the distribution of outward velocities 
towards a smooth form, as we increased the resolution of the 
simulations, led us to anticipate the converged form of this 
distribution. This was used as the basis for a more cohesive 
theoretical picture of the process, which takes into account 
both the local, momentum-conserving properties of the wind 
and the importance of the gravitational potential barrier 
which the winds must overcome. 

An example velocity distribution was then presented 
(Fig. [9| to demonstrate how such a modelling approach 
could work, and used to extrapolate predictions for the out- 
flow mass which were compared with the available simula- 



tion results (Fig. 10 1. The model was also compared with the 
observational estimates (Fig. 1 11 [ ) , demonstrating the poten- 
tial for the emergent global properties of these systems to 



12 This yield would be available from a supernova sweeping up 
a mean mass of 100 — 300Mq before a significant fraction of an 
initial energy of ~ 10 44 J was lost through cooling (depending on 
the IMF and structural factors of order unity) . 



be accounted for directly in terms of their internal structure 
and the energy sources which drive them. 
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APPENDIX A: SOFTENING LENGTH 



In section [3T4] there is a detailed comparison between simu- 
lation results generated using different numbers of particles. 
In order to make sure that this change in number is solely 
responsible for any difference in outcome, the gravitational 
softening length is crudely set to the same value for all cal- 
culations. 

The choice, lOpc, was made as it is appropriate for the 
higher resolution runs and these are results that are actually 
under the most scrutiny. However, this choice is inappropri- 
ate for the lower resolution runs. It is important to verify 
that the sensitivity of the results to particle mass that was 
exposed in Figs. [6] and [7] might not have been alleviated by 
more thoughtful choices. 

Some basic tests indicated that the results were indeed 
affected by changing the softening length, but confirmed 
that no choice could produce identical simulation results 
when run with very different particle mass. An example of 
one such test is shown in Fig. |A1| 
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Figure Al. Testing the consequences of choosing a more ap- 
propriate gravitational softening length for low-resolution simu- 
lations. This figure shows the time evolution (upper panel) and 
exit velocity distribution (lower panel) for a large disk galaxy at 
two contrasting resolutions, but with more appropriate choices of 
softening length than in Figs. [6] & [7] Whilst the low-resolutions 
results (solid lines) were indeed sensitive to this change, it clearly 
does not bring them in line with the high resolution calculations 
(dotted lines). In the upper panel, the colouring is as for Figs. [6] 
Red for ejected gas, magenta for the initial "wind", and blue 
for the ejected gas that has re-cooled back into the disk. 
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APPENDIX B: TIME STEPS 

The numerical methods used to illustrate this paper perform 
integration for each particle using time steps that are deter- 
mined by its local environment. As pointed out by Saitoh & 
Makino (2009), this will lead to very inaccurate integration 
in the case of strong shocks, where one locally-determined 
time step can be very long compared with the that of neigh- 
bouring regions. 

For the application to which the simulations have been 
put in this work, we are principally concerned with the ques- 
tion of whether this problem is responsible for the large 
change in results between low and high particle mass that 
is shown in Figs. [5] and [7] This possibility was ruled out 
by running the calculations again for the same initial con- 
ditions, but with the code modified such that the timestep 
for all particles is very short (about 5,000 years). 

The results of one such modified calculation is compared 
with the results of the standard scheme in Fig. |B1| The 
results are clearly affected by the use of much shorter time 
steps, but certainly not to the extent that they are changed 
by the use of smaller particle mass, as shown in Fig.s|6]and[7] 
Similar tests on other values of halo mass and particle mass 
produced similar outcomes, supporting the conclusion that 
rectifying the issues with time step allocation (as achieved 
by ?) would not completely remove the dependence of results 
on choice of particle mass. 
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Figure Bl. Testing the importance of timestep length in the 
low-resolution simulations. The dotted lines are the results of one 
of the resolution runs from Fig. [7] which used the standard cal- 
culation method where timesteps in the calculation are decided 
dynamically, based on the environment. The solid lines show the 
evolution of identical initial conditions where the calculation was 
changed only by enforcing the use of extremely short time steps 
throughout (Ai R2 5600 years). The key is as given, for definitions 
of the three components, see §3.1| 



